Quantum computers in phase space 
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We represent both the states and the evolution of a quantum computer in phase space using 
the discrete Wigner function. We study properties of the phase space representation of quan- 
tum algorithms: apart from analyzing important examples, such as the Fourier Transform 
and Grover's search, we examine the conditions for the existence of a direct correspondence 
between quantum and classical evolutions in phase space. Finally, we describe how to di- 
rectly measure the Wigner function in a given phase space point by means of a tomographic 
method that, itself, can be interpreted as a simple quantum algorithm. 
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I. INTRODUCTION 

Quantum mechanics can be formulated in phase space, 
the natural arena of classical physics. For this we can use 
the Wigner function [1], which is a distribution enabling 
us to represent quantum states and temporal evolution 
in the classical phase space scenario. In this paper we use 
the generalization of the familiar Wigner representation 
of quantum mechanics to the case of a system with a 
finite, A-dimensional, Hilbert space. Our main purpose 
is to develop and study the phase space representation of 
both the states and the evolution of a quantum computer 
(some steps in this direction were described in [2]). 

One can ask if there are potential advantages in using 
a phase space representation for a quantum computer. 
The use of this approach is quite widespread in various 
areas of physics (such as quantum optics, see [3] for a 
review) and has been fruitful, for example, in analyzing 
issues concerning the classical limit of quantum mechan- 
ics [4,5] . In answering the above question one should have 
in mind that a quantum algorithm can be simply thought 
of as a quantum map acting in a Hilbert space of finite di- 
mensionality (a quantum map should be simply thought 
of as a unitary operator that is applied successively to a 
system). Therefore, any algorithm is clearly amenable to 
a phase space representation. Whether this representa- 
tion will be useful or not will depend on properties of the 
algorithm. Specifically, algorithms become interesting in 
the large N limit (i.e. when operating on many qubits). 
For a quantum map this is the semiclas steal limit where 
regularities may arise in connection with its classical be- 
havior. Unraveling these regularities, when they exist, 
becomes an important issue which can be naturally ac- 
complished in a phase space representation. Therefore, 
having in mind these ideas, we conjecture that this rep- 
resentation may be useful to analyze some classes of al- 
gorithms. Moreover, the phase space approach may al- 
low one to establish contact between the vast literature 



on quantum maps (dealing with their construction, the 
study of their semiclassical properties, etc) and that of 
quantum algorithms. This, in turn, may provide hints to 
develop new algorithms and ideas for novel physics sim- 
ulations. As a first application of these ideas in this pa- 
per we examine several properties of quantum algorithms 
in phase space: We analyze under what circumstance it 
is possible to establish a direct classical analogue for a 
quantum algorithm (exhibiting interesting examples of 
this kind, such as the Fourier transform and other exam- 
ples which naturally arise in studies of quantum maps). 
We also shown that, quite surprisingly, Grover's search 
algorithm [6] can be represented in phase space and in- 
terpreted as a simple quantum map. 

To define Wigner functions for discrete systems various 
attempts can be found in the literature. Most notably, 
Wooters [7] proposes a definition that has all the desired 
properties (see below) only when A is a prime number. 
His phase space is an A x N grid (if A is prime) and a 
Cartesian product of spaces corresponding to prime fac- 
tors of A in the general case. Our approach here follows 
Wooters ideas and is closely related to that of Leonhardt 
[8] who analyzed Wigner functions for spin systems (both 
in the integer and half-integer case) and discussed in- 
teresting tomographic schemes to reconstruct the quan- 
tum state from measurements of marginal distributions. 
Other works connected to ours are those of Hannay and 
Berry [9] , who used discrete Wigner function in the con- 
text of studies of quantum chaos; Rivas and Ozorio [10] 
who define translation and reflection operators relating 
to the geometry of chords and centers on the phase space 
torus. Bouzouina and De Bievre [11] give a more ab- 
stract derivation of the same Wigner function related to 
geometric quantization. 

The paper is organized as follows: In Section II we 
briefly discuss the main properties of the Wigner repre- 
sentation in the continuous case. In Section III we review 
the main features of the discrete Wigner function. This 
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section is supposed to be self contained and to summarize 
known results (it also contains some original results and 
new explanations of old ideas). In Section IV we exam- 
ine the Wigner function of quantum states, some of which 
are of interest for quantum computation. In Section V we 
review the main properties of the temporal evolution in 
the Wigner representation. Here we describe some gen- 
eral results on the nature of temporal evolution in phase 
space and draw important analogies between quantum 
algorithms and maps. We explicitly analyze Grover's al- 
gorithm in phase space and also discuss the conditions 
for quantum evolution to have a direct correspondence 
with a classical map. In Section VI we consider the mea- 
surement of the Wigner function. We show that this can 
be done by means of a simple quantum computation that 
bears a remarkable similarity with a simulated scattering 
experiment. In Section VII we present our conclusions. 

II. CONTINUOUS WIGNER FUNCTIONS 

For a particle in one dimension, the Wigner function [1] 
is in one to one correspondence with the density matrix 
p and is defined as 

W(q,p) = J ^f Xvlh (<1 WW + A/2). (1) 

This function is the closest one can get to a phase space 
distribution for a quantum mechanical system. It's three 
defining properties are: (PI) W{q 1 p) is real valued, 
(P2) the inner product between states p\ and p 2 can 
be computed from the Wigner function as Tr[/5i/92] = 
2irh J dqdpWi(q,p)W2(q,p) and (P3) the integral along 
any line in phase space, defined by the equation aiq + 
o-iP = 0-3, is the probability density that a measurement 
of the observable a\Q + q 2 P has 03 as its result. This 
last property (the fact that W(q,p) yields the correct 
marginal distribution for any quadrature) is the most re- 
strictive one and, as Bertrand and Bertrand showed [12], 
together with P1-P2, uniquely determines the Wigner 
function. 

It is convenient to write W{q 1 p) as the expectation 
value of an operator, known as "phase space point op- 
erator" [7] (or Fano operators [13]). In fact, the Wigner 
function is 

W(q,p) = r Tv[pA(q,p)}. (2) 

The operators A(q,p) depend parametrically on q, p and 
can be written in terms of simetrized products of delta 
functions as: 

A(q,p) = :S(P-p)S(Q-q) : 

= [ dXdX ' c -ii(P-v)+i^(Q-ci) 
J (2^) 2 



where we have identified the continuous translation oper- 
ator £>(A, A') = cxp[--|(AP - X'Q)}. Thus, the above ex- 
pression shows that the phase space point operator is sim- 
ply the double Fourier transform of the phase space dis- 
placement operator (and, therefore, W(q,p) is the double 
Fourier transform of the expectation value of D(\, A')). 
It is even more convenient to rewrite these expressions 
as: 

A(q,p) = -LD(q,p)RD\q,p), (4) 

where R is the reflection (parity) operator that acts on 
position eigenstates as R\x) — \ —x). This means that the 
Wigner function is the expectation value of a displaced 
reflection operator. 

The proof of the three defining properties of the 
Wigner function (P1-P3) can be seen to follow from sim- 
ple properties of the phase space point operators. The 
fact that W(q,p) is real valued (PI) is a consequence of 
the hermiticity of A(q,p). Property (P2) follows from the 
completeness relation of A(q,p). In fact, one can show 
that these operators satisfy the relation 

Tr[A(q,p)A(q',p')} = -L_6(q - q')5(p - p'). (5) 

As a consequence of this, one can invert (2) and ex- 
press the density matrix as a linear combination of 
the phase space point operators. The Wigner func- 
tions determines the coefficients of such expansion: p = 
2nh J dqdpW(q,p)A(q,p). From this, the validity of the 
inner product rule can be easily demonstrated. The last 
property (P3) can be seen to be valid by observing that 
integrating A(q,p) along a line in phase space one always 
gets a projection operator. Thus, 

J dqdp5(aiq + a 2 p-a 3 )A(q,p) = \a 3 )(a 3 \, (6) 

where \a 3 ) is an eigenstate of the operator a±Q+a 2 P with 
eigenvalue a 3 . This identity can be shown by first writing 
the delta function as the integral of an exponential and 
then performing the phase space integration. We omit 
the proof here because we will present the corresponding 
one for the discrete case, which is done following very 
similar lines. 



III. DISCRETE WIGNER FUNCTIONS 

A. Preliminarities: discrete phase space 

We will consider here a quantum system with an TV 
dimensional Hilbert space (the case of a quantum com- 
puter is a specific example we will always have in mind 
but the formalism can be applied in other cases). In the 
Hilbert space we can introduce a basis B x = {\n),n — 
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0, . . . , N — 1} which we arbitrarily interpret as our dis- 
cretized position basis (with periodic boundary condi- 
tions: \n + N) = |n)). Clearly, the notion of a position 
observable may not have a direct physical interpretation 
in a discrete system but is a necessary ingredient if one 
wants to build a phase space representation. In the case 
of a quantum computer this basis can be, for example, 
the computational basis. Given the position basis B x , 
there is a natural way to introduce the conjugate mo- 
mentum basis B p = {\k), k = 0, . . . , N — 1} by means of 
the discrete Fourier transform. Thus, the states of B p 
can be obtained from those of B x as 

\k) = ^Y^^v(i2irnk/N)\n). (7) 

Therefore, as in the continuous case, position and mo- 
mentum are related by the Fourier transform. It is also 
important to recognize that the dimensionality of the 
Hilbert space of the system is related to an effective 
Planck constant. In fact, phase space should have a finite 
area (which we consider equal to one, in the appropriate 
units). In this area we should be able to accommodate N 
orthogonal states. If each of this states occupies a phase 
space area which is equal to 2nh, we see that: 

N=l/2wh. (8) 

In other words, N plays the role of the inverse of the 
Planck constant (and the large N limit is, in some way, 
the semiclassical limit). 

Once we have position and momentum basis we can 
construct the displacement operators in position and mo- 
mentum. Obviously, in this case it does not make sense 
to talk about infinitesimal translations. However, for dis- 
crete systems we can define finite translation operators 
U and V, which respectively generate finite translations 
in position and momentum [14]. The translation opera- 
tor U generates cyclic shifts in the position basis and is 
diagonal in momentum basis: 

U m \n) = \n + m), U rn \k) = exp{-2mmk / N)\k) , (9) 

(note that all additions inside kets are to be interpreted 
mod N). Similarly, the operator V is a shift in the mo- 
mentum basis and is diagonal in position: 

V m \k) =\k + m), V m \n) = exp(i2irmn/N)\n) (10) 

It is important to notice that the translation opera- 
tors U and V have commutation relations that directly 
generalize the ones corresponding to finite position and 
momentum translations in the continuous case. In fact, 
one can directly show the following identity: 

V P U 9 = U q V p exp{i2npq/N). (11) 

With these tools at hand we can introduce the dis- 
crete analog to the continuous phase space transla- 
tion operator D(q,p) = exp[— j^(qP — pQ)]- To do 



this we can first rewrite the phase space displacement 
as D(q,p) = cxTp(—iqP/h) cnp(ipQ/h) enp(ipq/2H) (an 
expression obtained by using the well known formula 
c A+B = e~[ A,B l/ 2 c A e B and the canonical commutation 
relation between Q and P). Therefore, by identifying 
the corresponding displacement operators, the discrete 
analogue of the phase space translation operator is: 

f(q,p) = U q V p exp(inqp/N). (12) 

These operators obey the simple composition rule 

f( qi ,pi)f(q 2 ,p 2 )=f(q 1 +q 2 ,p 1 + p 2 )e^ q - 

(13) 

where the phase appearing in the left hand side has a 
simple geometrical interpretation as the area of a triangle 
(see below). From the above definitions it is also simple 
to show that the phase space translation operators arc 
such that 

T*(q,p) = t(2N-q,2N-p) (14) 
= f(N-q,N-p)(-l) N+p+q . (15) 

Finally, it is worth mentioning the fact that the compo- 
sition rule can be used to show that 

f(Xq,Xp)=f x (q,p) (16) 

for any integer A, which is a very natural result. It may 
be tempting to define position and momentum operators 
also in the discrete case. In fact, we could attempt doing 
that by writing U and V as the exponentials of two opera- 
tors Q and P defined to be diagonal in B x and B p respec- 
tively: U = exp(-i2irP/N) and V = exp(i2TrQ / N) , with 

Q = J2 n n \ n )( n \ an d P = Sfcfc|fc)(^l- However, it turns 
out that these operators are not canonically conjugate 
since their commutator is not proportional to the iden- 
tity as is the continuous case (actually, it is well known 
that for finite dimensional spaces it is impossible to find 
two operators whose commutator is proportional to the 
identity). So, even though finite shifts in the discrete case 
have the same properties than their continuous counter- 
parts, this is not true for the position and momentum 
operators. 

It will also be useful to define a reflection (or parity) 
operator R as the one acting on states in the position 
(and momentum) basis as R\n) = \ — n) (again, the op- 
eration is to be understood mod N). It is simple to show 
that the parity operator obeys simple relations with the 
shifts U and V: 

UR = RU-\ VR = RV~ 1 . (17) 

It is important to recognize that the reflection operator 
is directly related to the Fourier transform. In fact, if we 
denote the operator implementing the discrete Fourier 
transform as Uft (i-e., the operator whose matrix ele- 
ments in the B x basis are (ti'|?7ft|i) = exp(i2imn' / N)) , 
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then the reflection operator is simply the square of the 
Fourier transform: 

R = U 2 FT . (18) 

Before ending this section it is worth mentioning other 
possibilities for the phase space description of finite quan- 
tum mechanics. The use of the discrete Fourier transform 
to relate position and momentum implies the imposition 
of periodic boundary conditions on both these variables, 
thus imposing the geometry of a torus on the phase space. 
This is by no means mandatory as the well known exam- 
ple of the angular momentum sphere -spanned by states 
\j, m),m = —j, ...j- illustrates. However, the ubiquitous 
appearance of the DFT in quantum algorithms somehow 
favors the choice of the torus as the preferred geome- 
try. It is also possible to slightly extend the scope of the 
present approach by imposing quasi-periodic boundary 
conditions, as is done in solid state chain models or in 
the treatment of the quantum Hall effect. This general- 
ization can be easily incorporated into the description of 
quantum algorithms [15] 

B. Discrete Wigner functions and phase space point 
operators 

To define a discrete Wigner function the most conve- 
nient strategy is to find the correct generalization of the 
phase space point operators A(q,p). Thus, we need to 
find a basis set of the space of hermitian operators (note 
that, as the Hilbert space is N dimensional such a ba- 
sis has iV 2 independent operators). We should do this 
in such a way that the resulting Wigner function gives 
the correct marginal distributions. However, generaliza- 
tion from the continuous to the discrete case has to be 
done with some caution. In fact, it is instructive to see 
how naive attempts to generalize expressions like (3) or 
(4) fail. In fact, generalizing equation (3) to the discrete 
case would lead us to define 

Ml,P) = ^2 E r(m,fc)exp{-i2 7 r (fcg ~ mp) }. (19) 

m,k— 

Unfortunately, this expression for A(q,p) turns out to be 
non-hermitian. This can be seen by taking the hermitian 
conjugate of the above equation and using the fact that 
the Hermitian conjugate of the phase space translation 
operator, as shown in (15), is such that T^(q,p) ^ T(N — 
q,N — p). Therefore, the above is not a good definition 
for the discrete phase space point operators. 

Another naive attempt to define discrete phase space 
operators is to generalize equation (4) which tells us that 
in the continuous case the operator A(q,p) is a displaced 
reflection. Thus, by writing 

Mq,p) = -^D{q,p)RD ] {q,p) - j^f(q,p)RI*(q,p), 

(20) 



one defines an operator that is hermitian by construction. 
However, in this case the problem is that the number 
of operators one defines in this way is not enough to 
build a complete set. In fact, if one uses the definition of 
T(q,p) given in (12) and then the commutation relations 
between U, V and R one finds: 

A(q,p) = ^U 2q RV- 2p exp{4Tripq/N). (21) 

As U and V are cyclic operators with period N, it is 
easy to see that if both q and p take values between 
and N the above expression only gives A 2 /4 independent 
operators. In fact, it is clear that the above equation 
implies that A(q + N/2,p) = A(q,p) and likewise with p 
so we only have N/2 x N/2 independent operators. 

Remarkably, the solution to the problems arising with 
the two above naive definitions is obviously the same: 
We should define, as it is done in the literature [9,11] 
the phase space point operators on a phase space grid 
of 2N x 2N points. This can be done simply replacing 
N — ► 2N in equation (19) or by taking half integer values 
for q and p in (21). The two definitions turn out to be 
equivalent. Therefore, the correct discrete phase space 
point operators are 

A(a) = j2Nf E r(A,A')exp{-i2 ^ q — P) } (22) 

= ^U q RV- p cxp(tn Pq /N) (23) 

where a = (q, p) denotes a point in the phase space grid 
with q and p taking values between and 2 A — 1. 

It is worth noticing that, as we have just defined phase 
space point operators on a lattice with 2N x 2N points, 
we have a total of 47V 2 such operators. However, it should 
be clear that these operators are not all independent. It 
is easy to verify that there are only N 2 independent phase 
space point operators since it can be proved that: 

A(q + a q N,p + a p N) = A(q,p) 

(24) 

for <J q ,o- p = 0,1. Therefore, it is clear that the A 2 phase 
space point operators corresponding to the first N x N 
subgrid of the phase space determine the rest. For the 
rest of the paper we will denote the first N x N subgrid 
as the set Gn (i-e., Gat = {a = (q,p);0 < q,p < N — 1}). 
The set G 2N will denote the full 2A x 2A grid. 

Before showing explicitly that the above operators en- 
able us to define a Wigner function with all the desired 
properties, it is worth pointing out some useful facts 
about A(a). In particular, it is convenient to realize 
that these operators are closely related to the phase space 
translation operators. In fact, one can show that by suc- 
cessive application of two phase space operators one al- 
ways gets a translation (this is also true in the continuous 
case) since 
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A(a)A(a>)=T(a-a>) — (25) 

Moreover, it is also useful to express the translation op- 
erator in terms of the A(a) operators by simply inverting 
the defition (22) and writing 



2JV-1 



f(n,k)= J2 A (l,P) exp(-i— (np-fcg)), (26) 

q,p={) 



a property that is also valid in the continuous case. 

From the above properties it is possible to show that 
A(a) are a complete set when a takes values in the first 
N x N subgrid Gn- Thus, taking the trace of (25) one 
gets: 

Tr[A(a)A(a')] = -L 5 N (q' - q)5 N (p' - p) (27) 

where both a and a' are in the grid Gn and Sn{q) = 
~k J2 n =o e 2mq / N is the periodic delta function (which is 
zero unless q = mod N). 

Therefore, according to all the previous arguments, the 
discrete Wigner function is defined as 



W(a) = Tr(i(a)p) 



(28) 



where a € G 2 n- Clearly, these 47V 2 values are not all 
independent since the Wigner function obeys the same 
relation than the phase space point operators: 

W(q + <7 q N,p + <J p N) = W(q,p) (-iyp1+°«P+°«°pN 

(29) 

As the operators A(a) are a complete set one can ex- 
pand the density matrix as a linear combination of such 
operators. It is clear that the Wigner function VF(a) are 
nothing but the coefficients of such expansion. Thus, one 
can show that 



p = 4N W(a)A(a) 

aEG N 

= JV ^ W(a)A(a). 



(30) 
(31) 



aeG 2 



The last expression can be obtained from (30) by notic- 
ing that the contribution arising from the four N x N 
subgrids are identical. 

Now it is simple to show that the Wigner function de- 
fined above obeys the three defining properties (P1-P3). 
The first one is a consequence of the fact that the op- 
erators A(a) are hermitian by construction. The second 
property (P2) is a consequence of the completeness of the 
set A(a) which enables one to show that 



Tt[ PiP2 }=N J2 W 1 (a)W 2 (a) 



(32) 



a€G 2 



Finally, the third property (P3) is more subtle and de- 
serves to be studied with some detail. The crucial point 
is to show that if we add the operators A(q,p) over all 
phase space points lying on a line L we always obtain 
a projection operator (this guarantees that adding the 
value of the Wigner function over all points in a line 
gives always a positive number, which can be interpreted 
as a probability). Before showing this, we should clearly 
define what do we mean by a line L{ni,n 2 ,ns) on our 
phase space grid. We will use the following definitions: 
A line L is a set of points of the grid defined as L — 
L(ni,n 2 ,n 3 ) = {(q,p) <G G 2 n, such that nip-n 2 q = n 3 , 
with < rii < 2N— 1}. Below, we will discuss some more 
details about the structure of lines on a grid G 2 m- Here, 
we only need to point out that one can also define a no- 
tion of parallelism between lines: two lines parametrized 
by the same set of integers n\ and n 2 are parallel to each 
other (in Figure (1) we show some examples of lines on 
a grid). 

So, let us show that by adding phase space point oper- 
ators on a line one always gets projection operators. We 
are interested in looking at the operator Al defined as 



Al= Mq,P)- 

(q,p)eL 



(33) 



It is clear that this operator can be rewritten as 

2N-1 



A L = y~] A(q,p) 5 2N {n\p-n 2 q-nz) 

q,p=0 

2N-1 2JV-1 



n 3 ) 



A=0 q,p=0 
2N-1 



^ J2 T A K,n 2 )e^A 



2N 



(34) 



A=0 



where to perform the sum over the phase space grid we 
used the fact that the Fourier transform of A(a) is given 
in (26). Now, as T(m,n 2 ) is unitary, it has N eigen- 
vectors \(f>j) with eigenvalues exp(—i2Tr<pj/N). Moreover, 
this operator is cyclic and satisfies T N = I. Therefore, 
as its eigenvalues are /Vth roots of unity, 4>j must be in- 
tegers. Then, we can rewrite equation (34) as: 



2JV-1 N 



^ = ^E E 



N 



A=0 j=0 



= ^£ 2 iv(2<Aj -n 3 ) \4>j){<t>j\- 

3=0 



(35) 



Therefore, it is clear that Al is a projection operator onto 
a subspace generated by a subset of the eigenvectors of 
the translation operator T(ni, n 2 ). The dimensionality d 
of this subspace is equal to the trace of Al- To calculate 
it, we just have to notice that in general 
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Tr[A(q,p)] 



1 

27V 



N-l 



i7T(q-2q')p/N 



(36) 



This last equation is easy to evaluate but the result de- 
pends on the parity of N. Thus, if TV is even then 
Tr (A(q, p)) = 1/N if both q and p are even and it is 
zero otherwise. On the other hand for odd values of N, 
one has that Tr(A(q,p)) = 1/2N for all values of q and 
p except when they are both odd where it has the oppo- 
site sign (i.e., it is equal to —1/2N). In the rest of the 
paper we will concentrate on analyzing the case of even 
N which is somewhat more relevant for a quantum com- 
puter (that has at least one qubit and therefore a Hilbcrt 
space which is even dimensional). Using the above result 
for the trace of A(a) one can see that the dimension- 
ality of the projector Al is simply given by 1/N times 
the number of points that belong to the line L which 
have even q and even p coordinates. An immediate con- 
sequence of this is that if n 3 is odd then d — (since the 
sum of two even numbers can never be odd). Finally we 
can write an explicit formula for d in the case of even N: 

1 N_1 

d =oE MAni)MAn 2 )e i3 33- (1 + (-l) n ») (37) 

1 A=0 

Let us illustrate the above result applying it to the sim- 
plest example: For a line L q defined as q — n 3 (i.e., ri\ = 
1, n i =0), the Wigner function summed over all points in 
L q ^ E M eL 9 W(Q,P) = E p W(n 3 ,p) = (n 3 /2\p\n 3 /2) 
if n 3 is even (and zero otherwise). Analogously, con- 
sidering horizontal lines (L p defined as p — n' 3 ) we get 
E {q , p) eL p W(q,p) = J2 q W(q,n' 3 ) = (n' 3 /2\p\n' 3 /2) if n' 3 
is even (and zero otherwise). These are just two examples 
a general result: this Wigner function always generates 
the correct marginal distributions (this, as in the contin- 
uous case, is the defining feature of W(q,p)). 

As a final point in this section where we reviewed 
some known (and other not so well known) results on 
discrete Wigner functions we think it is convenient to 
mention some simple properties of the phase space grid 
G2N ■ Most of the ideas we are using were introduced by 
Wootters [7]. As mentioned above we can simply define 
lines L{n\,ri2,nz) and introduce a notion of parallelism 
in the grid G2N- A foliation of the grid with a family of 
parallel lines is obtained by fixing ri\ and 712 and varying 
n 3 (in general, it is evident that two lines are parallel if 
the ratio n\/ni is the same). If N is a prime number 
then in the grid Gat (which has N x N points) there are 
exactly N(N + 1) distinct lines which can be grouped 
into N + 1 sets of parallel lines (these are N + 1 differ- 
ent foliations of the grid). If N is not prime or, as it is 
the case we are interested in, in the grid G2N this result 
is no longer true. For example, it is very clear that the 
equation mq — n 2 p = n 3 (mod 2N) has no solutions for 
odd values of n 3 and even values of n\ and ri2- So, it is 
not generally true that our lines have always 2iV points: 
As shown above sometimes they can have no points at all 



and in some cases one can construct lines with a number 
of points that is a multiple of 2N . This is the case if n\ 
and 712 have common prime factors with N . For example 
in the case 2N = 8, each line with m = ri2 = 1 have 8 
points, but lines with m = n 2 = 2 have no points if n 3 
is odd or 16 points each when n 3 is even. In any case, 
the simplest way to represent the lines in the phase space 
is by noticing that, as the topology of the lattice is that 
of a torus (due to the cyclic boundary conditions) lines 
wrap around the torus. The number of points in the line 
is related to the number of times it wraps around before 
it closes onto itself. Figure 1 shows an example of two 
sets of lines on the grid G2N for the case of N = 4. 

2q+2p=n 3 





q+p=n 3 
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FIG. 1. Left: lines satisfying the equation p + q = 713 (mod 
8) for values 713 = 0, 2. Each line has 8 points and there are a 
total of 8 distinct lines. Right: two of the four possible lines 
satisfying 2p+2q = 713 with 713 = 0, 2. In both plots horizontal 
(vertical) axis corresponds to position (momentum) basis. 

Let us summarize the results presented in this Section: 
We defined a Wigner function for systems with a finite 
dimensional Hilbcrt space (of arbitrary dimension N). 
The Wigner function is defined as the expectation value 
of the phase space operator A(a) given in (23). This 
definition is such that W(a) is real, it can be used to 
compute inner products between states and it gives all 
the correct marginal distributions when added over any 
line in the phase space, which is a grid G2N with 2A^ x 2A^ 
points. The size of the phase space grid is important to 
obtain a Wigner function with all the desired properties. 
The values of W(a) on the subgrid Gn are enough to re- 
construct the rest of the phase space (since the set A(a) 
is complete when a belong to the grid Gn)- However, 
the redundancy introduced by the doubling of the num- 
ber of sites in q and p is essential when one imposes the 
condition that all the marginal distributions should be 
obtained from the Wigner function. In the rest of the 
paper we will concentrate on applying this Wigner func- 
tion to study the states and the evolution of a typical 
quantum computer. 
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IV. WIGNER FUNCTIONS OF QUANTUM 
STATES 

To compute the Wigner function of any quantum state 
it is convenient to use equation (23) for A(q,p) and write: 

N-l 

W(q,p) = — n|p|n)e^^- g /2). (38) 

n=0 

Moreover, it is worth remembering that one only needs 
to compute iV 2 independent values (restricting to a = 
(q,p) € Gn) and from them one can reconstruct the re- 
maining 37V 2 ones by using (29). Before showing the 
specific form of the Wigner function of some states it is 
worth mentioning some general features of the Wigner 
function of pure states. In fact, if the quantum state is 
pure then p is a projection operator. Then, expanding 
p in terms of the phase space operators as in (31) and 
imposing the condition p 2 = p one gets 

W(a)=4N 2 ^ r(a,f3,j)W(f3)W(j). (39) 

where the function T(a,(3, 7), that depends on three 
phase space points (i.e, on a triangle) is 

r(a,/3, 7 ) = Tr(i(a)i(/?)i( 7 )) 

if either two or three of the points (a, (3, 7) have even q 
and even p coordinates. Otherwise T(a,(3, 7) = 0. In 
the above expression, which is valid for even values of 
N, S(a,(3, 7) is the area of the triangle formed by the 
three phase space points (measured in units of the ele- 
mentary triangle formed when the three points are one 
site apart). A similar expression is obtained by Wooters 
[7]. The three point function T(a, f3, 7) has then a simple 
geometric meaning and will play an interesting role in 
determining the properties of the temporal evolution (as 
Unitary evolution preserves pure states, it will be rep- 
resented by a linear map in phase space that leaves T 
invariant, as shown see below). 

A. Position and momentum eigenstates, and their 
superpositions 

Here, we first evaluate the Wigner function of a po- 
sition eigenstate (a computational state of the quantum 
computer) p qo = \qo)(qo\- It is straightforward to obtain 
a closed expression for W(q,p): 

W go (q,p) = l^( qo \UiRV-P\ qo }e^/ N 

= ±-6 N (q-2q ) (-l)K«-2 9 o) N (41) 



where zn denotes z modulo N. This function is zero 
only on two vertical strips located at q = 2q modulo N. 
When q = 2q , W(q,p) takes the constant value 1/2N 
while for q = 2q ± N it is 1/2N for even values of p 
and —1/2N for odd values. These oscillations are typical 
of interference fringes and can be interpreted as arising 
from interference between the q = 2q strip and a mir- 
ror image formed at a distance 2N from 2qo, which is 
induced by the periodic boundary conditions. The fact 
the the Wigner function becomes negative in this inter- 
ference strip is essential to recover the correct marginal 
distributions. Adding the values of W(q,p) along a ver- 
tical line gives the probability of measuring q/2 which 
should be 1 for q = 2qo and zero otherwise. A very sim- 
ilar calculation can be done for a momentum eigenstate 
p = \k )(k \. The result is very similar except that now 
the strips are horizontal. 

It is interesting to analyze also the Wigner func- 
tion of a state which is a linear superposition \ip) = 
{\Qo) + e ~ l( ^\qi}) /V2. Again one can simply get a closed 
expression for W(q,p) which is 

W(q,p) = \ [W qo (q,p) + W qi (q,p) + AW qo , qi (q,p)\ 

where the interference term is 

AW qo , qi (q,p) ^^S N (q) (-1)* cos (^P + ^ , 

with q = (jo + qi — q and A = 2N/(qo — qi). Thus, this 
Wigner function has two direct terms which simply cor- 
respond to the two computational states and an interfer- 
ence term which is peaked on the vertical strips located 
at the midpoint q = qa + qi mod N. On this strip the 
Wigner function oscillates with a wavelength that is in- 
versely proportional to the separation between the main 
fringes. This oscillatory pattern has its corresponding 
mirror image and the interference between them produce 
some more oscillations. In figure 2 we show the Wigner 
function for a state and a superposition of two such com- 
putational states. The plot is much more eloquent than 
any equation and shows both the presence of the main 
peaks, its mirror images and (in the case of a superposi- 
tion state) the interference fringes. 
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FIG. 2. Wigner function of a position eigenstate (left) and 
a linear superposition of two position eigenstates (right). Hor- 
izontal (vertical) axis correspond to position (momentum) ba- 
sis. For the computational state |qo), the Wigner function is 
positive (black) on a vertical strip located at q — 2qo and it 
has negative (white) values on another vertical strip located 
at q = 2qo ± N. For a superposition of two computational 
states (right), the Wigner function also oscillates on the strip 
located at q = go + qi with a wavelength that depends upon 
the distance between the two interfering states. 



B. Other quantum states 

The Wigner function can be computed in closed form 
for some other interesting states. In figure 3 we display 
such function for the completely mixed state (that, as 
we mentioned above, is nonzero only when both q and p 
are even). We also show the Wigner function for a pure 
state constructed as a Gaussian superposition of compu- 
tational states (with periodic boundary conditions, i.e. a 
sum of a Gaussian state centered about a phase space 
point and all its infinite mirror images). The Gaussian 
wavepacket has a Wigner function with a positive Gaus- 
sian peak and shows three other Gaussian peaks that are 
modulated by interference fringes. They correspond to 
the interference between the main peak and its mirror 
images (notice the different orientation of these fringes). 



p(t + l) = U P (t)U f . 



(42) 



FIG. 3. Left: Wigner function for the completely mixed 
state. It is positive (black) only in the phase space points with 
even values of q and p being zero (grey) elsewhere. Right: 
Wigner function for a periodic Gaussian wavepacket. This 
function has a positive (black) peak and three regions where 
it displays oscillatory behavior. These regions can be under- 
stood as being originated from the interference between the 
positive Gaussian peak and its mirror images created by the 
periodic boundary conditions. Horizontal (vertical) axis cor- 
respond to position (momentum) basis. 



V. TEMPORAL EVOLUTION IN PHASE SPACE 

Temporal evolution of a quantum system can also be 
represented in phase space. If U is the unitary evolution 
operator that takes the state of the system from time t 
to time t + 1, the density matrix evolves as 



Using this, it is simple to show that the Wigner function 
evolves as 

W(a,t + 1)= J2 Z aP W(j3,t). (43) 

where the matrix Z a p is defined as 

Z aP = N Tr{A{a)UA{(3)U r ). (44) 

Therefore, temporal evolution in phase space is rep- 
resented by a simple linear transformation (which is, 
of course, an immediate consequence of linearity of 
Schrodinger equation). Unitarity imposes some con- 
straints on the matrix Z a p: As purity of states is pre- 
served, temporal evolution must preserve the structure 
of the constraint equation (39). Therefore, the matrix 
must leave invariant the three point function T(a,(3, 7), 
i.e. 



r(a',/3', 7 ') 



E 



Z a 'aZf3if}Zy j r(a, (3, 7). 



(45) 



The real matrix Zap contains all the information about 
the nature of temporal evolution. In general, this matrix 
connects a point a with many other points (3. Therefore, 
temporal evolution will be generally nonlocal in phase 
space. This is a unique quantum mechanical feature: In 
fact, classical systems evolve in phase space following a 
flow of classical trajectories. In such case, the value of 
the classical distribution function W(a, t+l) is equal to 
the value W(/3, t) for some point (3 which is a well defined 
function of a and t. One may ask what kind of Unitary 
operators generate a local dynamical evolution in phase 
space. Below, we will give examples of this type. 



A. Phase space translations 

It is simple to show that if we consider the unitary 
evolution to be the phase space translation operator 
T(a) = T(q,p) then temporal evolution is local in phase 
space. Moreover, quantum and classical evolution are 
identical in this case since the value of the Wigner func- 
tion is rigidly translated in phase space: 



U = T(a) 



W{a,t + 1) = W{a-2a,t). (46) 



Notice that the factor of 2 in the above equation is simply 
a consequence of the fact that we are working in a phase 
space grid of 2N x 2N. Another example of a local evolu- 
tion is the one associated with the phase space operators 
A(a) themselves. In such case the resulting evolution is 
not just a translation but a translation combined with a 
reflection: 



U = A(a) 



W(a,t + l) = W(2a-a,t). (47) 
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The above are two simple examples of families of unitary 
operators with corresponding local phase space represen- 
tations. As we will see below, the family of operators 
with a direct correspondence between quantum and clas- 
sical evolution has other interesting members. 

B. Fourier transform 

The discrete Fourier transform is a unitary operator 
that is widely used in quantum algorithms [6]. As posi- 
tion and momentum basis are related by this operator one 
expects that it should have a rather simple phase space 
representation. In fact, this is indeed the case: One can 
easily show that the Fourier transform is represented as 
a rotation in 90 degrees: 

U = U FT W(q,p, t + 1) = W(-p, q, t). (48) 

where — p is the additive inverse mod(iV). Thus, the 
Fourier transform is also represented by a local opera- 
tion in phase space (and, in this sense it is a completely 
"classical" operation, acting independently on each of the 
N x N subgrid). For example, by applying the Fourier 
transform to the quantum states whose Wigner functions 
are shown in Figure 2 (a computational state and a su- 
perposition of two such states) one obtains the resulting 
Wigner function by rotating the above figure by 90 de- 
grees (i.e., one gets a momentum eigenstate or a super- 
position of two such states where the vertical pattern in 
Figure 2 is mapped into the same horizontal pattern). 
It is clear that by applying the Fourier transform twice 
one gets a rotation by 180 degrees. This is nothing but 
a reflection, which is what equation (18) tells us. 

C. Quadratic (cat) maps propagate classically 

The above operators (rigid translations, reflections and 
the Fourier transform) are just some examples of unitary 
operators which generate a classical evolution in phase 
space. A more general family with such property will be 
discussed here. This family is formed by the quantization 
of classical dynamical systems with a quadratic Hamil- 
tonian and a finite phase space with periodic bound- 
ary conditions (as the phase space is a torus, the classi- 
cal transformations are the linear automorphisms of the 
torus [16]). We will briefly present these unitary opera- 
tors here and show that classical and quantum evolution 
coincides (a fact that has been shown before, using differ- 
ent techniques, by Hannay and Berry [9]). Consider the 
following two-parameter family of operators (it is possible 
to consider a slightly more general case, with one more 
parameter, but to simplify the presentation we restrict 
to this simpler case [15]). 

U cat =V b TV a . (49) 



The operators V a and T are respectively diagonal in the 
position and momentum basis and satisfy 

Vj\n) = cxp(-i27rn 2 (l - j)/2N)\n) 
T\k) = cxp(-i2 7 rfc 2 /2A r )|fc) (50) 

where j is an integer. The operator U cat can be sim- 
ply thought of as the evolution operator of a kicked sys- 
tem with a Hamiltonian in which kinetic and momentum 
terms are turned on and off alternatively. The parame- 
ters a and b are integers that measure the strength of the 
potential kicks. It is straightforward to find out what is 
the classical system to which the above operator corre- 
sponds: One way to do this is to take the matrix elements 
of (49) in the computational basis and show that 

(n'\U cat \n) = K cxp(i27r(an 2 + bn' 2 - 2nn')/2N) (51) 

where if is a normalization constant. The exponent in 
this equation can be interpreted as the classical action of 
the system, where n and n' the final and initial values 
of the coordinates. Therefore, the classical equations of 
motion corresponding to this system are 

n = bn'+p' p = (ab - l)n' + ap' . (52) 

This classical system has been extensively studied [16]: 
For integer values of a and b, it is a member of the famous 
family of "cat" maps [16], i.e. all linear automorphisms of 
the torus. The system is chaotic when the eigenvalues of 
the linear transformation M mapping a' = (n' 7 p') onto 
a = {n,p) as in (52) are both real (otherwise the map is 
integrable). This is the case when when TrM. = a+b > 2 
(notice that M. has unit determinant). In particular, 
when a — 2 and 6=1 this is the so-called Arnold- 
cat map (n = n' + p', p = n' + 2p. This special cat, 
according to (49), simply has a kinetic kick followed by 
a potential kick where the potential is an upside-down 
harmonic oscillator. 

The reason why quantum and classical evolution are 
identical in this case can be shown by using our previous 
results. In fact, we need to compute the matrix Z a p that 
evolves the Wigner function. To do this calculation it is 
useful to first notice that, if the evolution operator U is 
the one given in (49), then the following identity holds: 

U cat A(a) = A(Ma) U caU (53) 

where, as above, the linear transformation is the one 
given in (52). This means that the unitary evolution op- 
erator maps the phase space point operator in the same 
way as the classical dynamics does with the phase space 
point. Using this, it is simple to show that the Wigner 
function evolves classically: 

U = U cat W(a,t+l) = W(M- 1 a,t). (54) 

Our result above is not so surprising when seen from 
the perspective of ordinary quantum mechanics of con- 
tinuous systems. In fact, it is well known that quantum 
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and classical evolutions are identical in phase space if 
and only if the classical system has a Hamiltonian that 
is quadratic in p and q (i.e. a harmonic oscillator). Here, 
we showed that the same result is valid for systems with 
a finite dimensional Hilbert space. It is worth pointing 
out that it is possible to design a simple quantum circuit 
to implement the evolution operator U cat [17]. This can 
be done by noticing that the potential kick is diagonal in 
the computational basis and to implement it one simply 
needs a circuit with the same complexity than the Fourier 
transform (in this way one can construct a circuit imple- 
menting a controlled phase gate where state \n) acquires 
a phase that depends quadratically in n [15]). The opera- 
tor corresponding to the kinetic kick con be implemented 
by a similar circuit in between a Fourier transform and 
its inverse. 



D. Boolean gates in phase space 

The family of quantum cat maps is rather large but 
does not contain some operators that are more natural 
from the point of view of quantum computation. We 
study here such operators by analyzing when does quan- 
tum and classical evolution are identical in phase space if 
the dynamics is generated by an arbitrary boolean gate. 
This kind of operation, which is a permutation of the 
computational basis, is an essential ingredient in Shor's 
factoring algorithms and others [6]. So, let us consider a 
one to one function / : Zn — * Zn (/ is a permutation 
of the first TV integers). Given / we can define a unitary 
operator Uf whose action in the computational basis is 
to permute vectors in the same way / permutes numbers: 

U f \n) = \f(n)) (55) 

The Boolean function / corresponds to some general 
(reversible) classical circuit. Moreover, one can also as- 
sociate with / a classical map in the phase space. Such 
map is the one permuting the vertical strips correspond- 
ing to the different positions according to the function 
/, i.e. it maps the vertical strip labeled by n into the 
one labeled by / (n)). Here, we want to determine under 
what condition the phase space evolution associated with 
Uf is identical to the classical one (i.e., to the mapping 
of the vertical strips). The answer to this question is 
remarkably simple: the quantum map is identical to the 
classical one if and only if f(n) = n + a (mod(iV), for any 
integer a). That is to say, quantum and classical permu- 
tations are identical only in the case of / corresponding 
to a cyclic shift. 

The proof of the above statement is simple: We will 
consider the action of Uf on two classes of states. First 
we will look at how U / transforms the Wigncr func- 
tion of a computational state (i.e., a density matrix like 
p = \no)(no\, shown in Figure 2 (a)). Then, we will 
analyze how U / modifies the Wigner function of the in- 
terference term between two computational states (i.e., a 



density matrix of the form p = |no)(ni| ± |«i)(no| shown 
in the interference pattern seen in Figure 2b). These 
operators form a complete basis of the space of hermi- 
tian matrices. Therefore, by looking at the evolution 
of each of these states we are also able to understand 
the evolution of the most general state. As Uf simply 
permutes computational states, it is straightforward to 
show that the Wigner function of any such state (shown 
in Figure 2a) will be transformed according to the clas- 
sical map, i.e. the resulting Wigner function is obtained 
from the original one just by permuting vertical strips 
according to the classical map 2n — > f{2n)). However, 
the fate of a state corresponding to the superposition be- 
tween two computational states is drastically different: 
As shown above, the interference fringes are originally lo- 
cated at the midpoint between the two interfering states 
(i.e. at point n = n + n\). Therefore, for the state to 
evolve classically one needs to impose that this midpoint 
is mapped into the midpoint of the transformed states 
(i.e. that 2/(m + n 2 ) = f(2n 1 ) + f(2n 2 )). Moreover, 
the original state has oscillations with a wavelength that 
depends on the distance between the two computational 
states (A = 27V/ (no — n\). Therefore, for the state to 
evolve classically these fringes must remain unchanged 
and the only possibility for this to happen is that the 
distance between the states remains constant (i.e. that 
2(n — Tii) = f(2n ) — /(2m)). The only solution for 
these two constraints is that f(n) = n + a, i.e. a linear 
function corresponding to a rigid shift in the computa- 
tional basis. 

This result can be generalized as follows: Instead of 
defining Uf as in (55) we can consider a more general 
operator Uf. g which is now associated with two inte- 
ger one to one functions in such a way that Uf. g \n) — 
exp(i2ng(n)/N)\f(n)). Following the same steps de- 
scribed above one can show that the operator acts in the 
same way as the corresponding classical map if and only 
if both functions are such that f(n) = n + a, g(n) = n + b. 
In such case this operator is nothing but the phase space 
translation T(a,b). Thus, our result shows that if one 
defines the quantum analog of a classical gate as we did, 
the only such unitary operator which generates the same 
phase space dynamical evolution as its classical analog is 
the one corresponding to a linear rigid shift. 

It is interesting to notice that the structure of the 
above proof (that focuses on how quantum interference 
terms are affected by quantum or classical evolution) 
could also be applied to ordinary quantum mechanics (of 
continuous systems) and sheds new insight on the reason 
why only linear systems evolves classically for all times. 
We remark that this result has remained unnoticed for 
quantum computers and may offer some help in our ef- 
forts to understand the differences between these systems 
and their classical counterparts. 



10 



E. Nonlocal (quantum) evolution in phase space 

Apart from the above simple examples with direct clas- 
sical counterpart, most of unitary operators generate a 
rather complicated nonlocal evolution for the Wigncr 
function. As an example we can explicitly compute ma- 
trix elements of Z a p for some simple case: If we consider 
a bit flip of the least significant qubit (which is associated 
with cri°^) we can show that: 

U = (4 0) Z a0 = Z 0a = 2/N if a is even 

(56) 

and is equal to zero otherwise (a = (q,p) is even when 
both q and p are even). Other matrix elements of Z a p 
are simple to compute but this is enough for the purpose 
of exhibiting the nonlocal nature of the evolution. The 
above equation implies that when U = <r^ the next value 
of W(0) is proportional to the sum of V^(a) over all even 
values of a which is indeed a rather nonlocal mapping. 

Another simple example of nonlocal evolution is con- 
structed by considering a unitary operator U^J 2 " 1 which 
consists of doing nothing to the most significant qubit 
and applying the Fourier transform to the rest (this is 
a Fourier transform on a space of dimension N/2). In 
such case, the resulting evolution bears some relation to 
a 90 degrees rotation (as the complete Fourier transform 
does) but is a rather nonlocal beast. It can be seen that 
this operation can be understood as the quantization of 
a Bernouilli shift. In fact it is an essential ingredient in 
the construction of the quantum bakers map (see [18-20] 
for more details). Here we are only concerned with the 
phase space description of this operation. Classically, 
when operating on a vertical strip, i.e. a computational 
state - this circuit shifts it by n' — > 2n(modN), contracts 
it to half height (and double thickness) and rotates by 
90°. Remnants of this behavior can be seen on the left 
plot of Figure 4 where the operation has acted on the 
computational state illustrated in Figure 2. Notice that 
the phase space representation shows strong diffraction 
effects which lead to a non local Z a p. It is not diffi- 
cult to understand where these effects come from: the 
unchanged qubit divides the phase space in two regions 
and the map acts discontinuously on them, thus creating 
the same effect as a screen in an optical system. This is 
contrasted in the right part of Figure 4 with the action 
of Arnold's cat map on the same state. In such case, the 
vertical strip is simply transported by the classical linear 
flow. 
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FIG. 4. Diffraction effects in the Wigner propagation. Hor- 
izontal (vertical) axis correspond to position (momentum) ba- 
sis. In the left plot we show the resulting Wigner function ob- 
tained by applying the operator U^/ (corresponding to a 
Bernoulli shift) to a computational state. The classical prop- 
agation would contract and rotate the line shown in Fig.2, but 
the quantum propagation shows strong diffraction effects. In 
the right we show the result of propagating the same initial 
state with the 'cat' map (equation (49)). In this case the 
Wigner function evolves classically and shows no diffraction 
effects. 



F. Quantum algorithms in phase space 

Quantum algorithms are nothing but unitary operators 
that have been designed so that after applying them sev- 
eral times, they produce a state that encodes the answer 
to some computational problem. A measurement on such 
state gives us information about this answer. Quantum 
algorithms are designed in such a way that they have a 
regular behavior as function of N. Thus, at least in some 
cases, one should be able to observe this regularities in a 
phase space representation. In fact, such representation 
becomes more useful also in the large N limit, which is 
the semiclassical limit (since 1/N is an effective Planck 
constant). The study of the phase space representation 
of quantum algorithms is only beginning but there are 
already some interesting results showing that some al- 
gorithms may find phase space as a natural arena to be 
represented in. The example we would like to point out 
here is Grover's search algorithm whose phase space rep- 
resentation is shown in Figure 5 
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FIG. 5. Grover's algorithm in phase space. Horizontal (ver- 
tical) axis correspond to position (momentum) basis. After 
only five iterations the initial state (a momentum eigenstate) 
is transformed into a position eigenstate located precisely at 
the position of the marked item (see text). 

In this Figure we display the Wigner function of the 
quantum state of a computer after every iteration of 
Grover's search algorithm. Our system has an N = 32- 
dimensional Hilbert space (i.e., the computer has just 5 
qubits) and the algorithm is designed to search for the 
marked item which we chose here to be q — 16 (indicated 
with an arrow in the plot). The initial state is chosen to 
be an equally weighted superposition of all computational 
states which, for the purpose of making the plots more 
visible, we chose to be a non-zero momentum state (the 
usual choice for the initial state in Grover's algorithm 
is a zero momentum state but the algorithm works as 
well with an initial state |fc) with nonzero momentum). 
The quantum search algorithm consists of the iteration 
of a unitary operator Ug which can be decomposed as 
Ug = UkU - The operator U represents a call to an or- 
acle that only distinguishes the marked item ui from the 
rest: U Q = I — 2\oj)(lo\. The operator Uk (the inversion 



about the mean) has no information about the marked 
item and is chosen as Uk = I — 2|fc)(fc|. The quantum 
computer should be initially prepared in the state Af- 
ter a number of iterations T tt^/N /4 the probability of 
detecting the position q = lu becomes of order unity (with 
small correction decaying as 1/N). Thus, by measuring 
the position at this time one would discover the marked 
item. From the figure one clearly observes that this al- 
gorithm is very simple when seen from a phase space 
representation. The initial state has a Wigner function 
which is a horizontal strip (with its oscillatory compan- 
ion). After each iteration W(q,p) shows a very simple 
Fourier-like pattern and becomes a pure coordinate state 
at the end of the search (in our case, the optimal number 
of iterations is T w tt^/N/4: w 5). This representation 
shows that, as a map, Grover's algorithm has a fixed 
phase space point with coordinate equal to the marked 
item (q — 16 in our case) and momentum equal to the 
one of the initial state. 

It is interesting to speculate if this phase space repre- 
sentation may be of guidance to create new quantum al- 
gorithms or improve on existing ones. We cannot present 
clear evidence to support this but we believe this will be 
the case. However, some careful thought is required. At 
this moment we can only offer the reader some of our 
experience on the way one should not proceed: View- 
ing Figure 5 one may be tempted to try to find a quan- 
tum map that has the same fixed point as the above 
but transform horizontal strips into vertical ones faster 
than Grover's. It is interesting to notice that such maps 
do exist. In fact, one just needs a hyperbolic (instead 
of an intcgrablc) map with the same fixed point struc- 
ture. In fact, a hyperbolic map with the same fixed point 
structure would perform the same operation exponen- 
tially faster, because it would transform a state initially 
along the stable direction into one along the unstable one 
in a time of order log(iV) . Such hyperbolic maps can be 
easily constructed as variations of the quantum bakers 
map. However, their construction requires from the or- 
acle much more information than a simple yes/no query 
and thus do not qualify as fast search algorithms. 

VI. MEASURING THE WIGNER FUNCTION 

In this section we review a general procedure to di- 
rectly measure the value of the Wigner function of a 
system in any phase space point (the procedure applies 
for discrete and also for continuous systems). This mea- 
surement strategy was originally proposed in [21] as an 
application of a very general tomographic scheme. It 
generalizes previous proposals that were put forward to 
measure the Wigner function in quantum optics or cavity 
QED setups (see [3] for a review). 

The basic idea of the measurement is best described 
in terms of the quantum algorithm represented by the 
simple quantum circuit shown in Figure 6. It works 
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as follows: Suppose you want to measure the Wigncr 
function of a system that is prepared in some quantum 
state p. What you need to do is to bring this system in 
contact with an ancillary qbit prepared in the state |0) 
and then "run" the quantum computation described by 
the circuit shown in Figure 6. The ancillary qbit plays 
the role of a "probe particle" in a scattering-like experi- 
ment. The algorithm represented in the circuit is: i) Ap- 
ply an Hadamard transform to the ancillary qbit (where 
H\0) = (|0> + H\l) = (|0) - |1))/V2), ii) Ap- 

ply a "controlled-C/" operator (if the ancilla is in state |0) 
this operator acts as the identity for the system but if the 
state of the ancilla is |1) it acts as the unitary operator U 
on the system), iii) Apply another Hadamard gate to the 
ancilla and finally perform a weak measurement on this 
qbit detecting its polarization (i.e., measuring the expec- 
tation values of Pauli operators a z and cr x ). It is easy to 
show that the above circuit has the following remarkable 
property: 



Re{Tr{Up)), 



-Im{Tr(Up)). (57) 



Thus, the final polarization measurement reveals a prop- 
erty determined both by the initial state p and the uni- 
tary operator U. 
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FIG. 6. Circuit for measuring Re(Tr(U p)), for a given uni- 
tary operator U. 

It is clear that this circuit can be used for dual pur- 
poses: On the one hand, we can use it to extract infor- 
mation on the operator U if we know the state p. On the 
other hand, the same circuit can be used to learn about 
the state p by using some specific operators for U. In 
this sense, this circuit can be adapted to act either as a 
tomographer or as an spectrometer which, therefore, are 
two dual faces of the simple quantum computation shown 
above [21]. The analogy with a typical scattering exper- 
iment is quite clear: the probe particle is prepared in a 
well determined state, then it interacts with the system. 
From the measurement of some statistical distribution of 
the probe (in this case the polarization) we are able to 
infer properties of either the interaction or the state of 
the scatterer. 

In [21] we showed that this circuits can be used to 
measure the Wigner function. This is done by simply 
choosing U = A(a). By doing this, every time one runs 
the algorithm a direct measurement of W(a) is obtained 
(for any given phase space point a). This measurement 
strategy can be, in principle, used both for the contin- 



uous and the discrete case. For example, in the contin- 
uous case, to measure W(a) of a quantum particle one 
has to run the algorithm with the system in the initial 
state D(a)pD^ (a) (obtained by displacing p) and using 
U = R (it is clear that this is equivalent, but simpler, 
than applying the circuit directly with U — A(a)). The 
measurement of the Wigner function for the discrete case 
can be done in a similar way but in that case it is more 
convenient to directly use U = A(a). In both cases, as 
the operator A(a) is both unitary and hermitian, the 
measurement of the Wigner function only requires deter- 
mining the z-component of the polarization. 

Measuring directly the Wigner function of a quantum 
system has been the goal of a series of experiments pro- 
posed and performed in various areas of physics (all deal- 
ing with continuous systems, [22]). It is interesting to 
notice that the above measurement strategy precisely de- 
scribes the experiment recently performed to determine 
the Wigner function of the electromagnetic field in a cav- 
ity QED setup [23,24]. In fact, in such case the system is 
the mode of the field stored in a high-Q cavity and the 
ancillary qbit is a two level atom. Remarkably, this ex- 
periment can be easily described in terms of the applica- 
tion of the algorithm shown in Figure 6: i) The atom goes 
through a Ramsey zone that has the effect of implement- 
ing an Hadammard transform. A radio-frequency source 
is connected to the cavity displacing the field in phase 
space (by an amount parametrized by a), ii) The atom 
goes through the cavity interacting dispersively with the 
field. The interaction is tuned in such a way that, if 
the atom is in state \g) nothing happens but if the state 
of the atom is |e) the state acquires a phase shift of it 
per photon in the cavity. So, this interaction is simply a 
controlled-exp(— iwN) gate (where N is the photon num- 
ber operator) which is nothing but a controlled reflection, 
iii) The atom leaves the cavity entering a new Ramsey 
zone that enforces a new H-gate. Finally the atom is 
detected in a counter either in the ground \g) or in the 
excited |e) state. The experiment is repeated many times 
and the Wigner function is obtained as the difference be- 
tween the probability of the atom being in the excited 
and ground states: W(q,p) = 2(P(e) — P(g))/h. As we 
see, this cavity-QED experiment is an example of a tomo- 
graphic experiment represented by the circuit of Figure 
6. The important point is that the method is generaliz- 
able to arbitrary systems (continuous or discrete) when 
viewed in terms of the above circuit. 

For the discrete case, one can show that the circuit in 
Figure 6 can be efficiently decomposed as a sequence of 
elementary gates. For this, one only needs to implement 
controlled-?/, V and R operations. All these can be im- 
plemented via efficient networks like the ones shown in 
[25] which require a number of elementary gates that de- 
pends polinomially on log(N). For example, to decom- 
pose V as a sequence of one qbit gates we can notice that 
its action on a computational state \n) is: 
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V p \n) = Yl [exp(i27rp2V7V)]™ , \n) 



(58) 



i=0 



where is the binary expansion of the number n. From 
this expression it is clear that to implement V we just 
have to act independently on each qbit. The operator 
(pi acting on the i-th qbit should be such that <j>i\0) = 
|0) and <f>i\l) = cxp(i2np2 l /N)\l). On the other hand, 
to implement U, which is a shift in the position basis, 
we can either use a circuit which performs addition mod 
N as described in [25] or notice that a circuit for U is 
obtained from the one corresponding to V by using U — 
Up^VUft- Finally, to implement the reflection operator 
we use the identity (18) which states that two Fourier 
transforms are equivalent to a reflection. Therefore, the 
complete network can be obtained straightforwardly. 

For small values of N the resulting networks are indeed 
very simple. In [21] we have presented the results of the 
measurement of the Wigner function for the particular 
case of N = 4 (two qbits) for a variety of initial states. In 
this case the reflection operator is simply a controlled not 
(CNOT) gate (where the control is in the least significant 
qbit). The circuit for U is the same CNOT followed by 
bit flip in the control. Analogously, the circuit for V is a 
sequence of controlled phase gates. Thus, the complete 
measurement circuit has at most one Tofolli gate and 
several two qbit gates. We experimentally measured the 
Wigner function of the four computational states of a 
two qbit quantum computer. The ideal result is shown 
in Figure 2 and the measured Wigner functions are seen 
in 7. 

To perform these experiments we used a liquid sam- 
ple of trichloroethylene (TCE) dissolved in chloroform. 
This molecule has been used in several three qubit exper- 
iments where the proton (H) and two strongly coupled 
13 C nuclei (C\ and C 2 ) store the three qbits [26]. In our 
case, we used C\ as our control qbit (the probe particle 
in Figure 6) and the pair H-C2 to store the state whose 
Wigner function we tomographically measure. The mea- 
sured coupling constants are JhCi = 200.76 Hz, Jhc 2 = 
9.12 Hz and JciC 2 = 103.06 Hz while the C1-C2 chemical 
shift is Sc 1 c 2 — 908.88 Hz. To evaluate the Wigner func- 
tion in each of the independent 16 phase space points we 
expressed the algorithm of Figure 6 as a sequence of r.f. 
pulses and delays (the number of pulses in each sequence 
depends on the phase space points and varies between 
5-20 and take at most 100ms to execute (<C T\,T% of 
our sample)). Experiments were done at room tempera- 
ture and we used temporal averaging to distill the initial 
pseudo-pure states [27] . The experiments were done on a 
standard 500Mhz NMR spectrometer (Bruker AM-500 at 
LANAIS in Buenos Aires and DRX-500 at Los Alamos). 
We used a 5mm probe tuned to 13 C and 1 H frequencies of 
125.77 MHz and 500.13 Mhz, respectively. The measured 
Wigner functions shown in figure 7 agree with the ideal 
one within (in the worst case) approximately 15% error. 
The most important sources of errors are understood and 



come from the effects of strong coupling (that alter the 
behavior of our gates) and numerical uncertainty in in- 
tegrating the spectra. It is clear that the result of this 
simple experiment, that illustrate the tomographic mea- 
surement of a discrete Wigner function, agrees very well 
with the theoretical expectation. 
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FIG. 7. Measured Wigner functions for the four computa- 
tional states of a two-qbit system (built with a liquid sam- 
ple of TCE in an NMR spectrometer). Horizontal (vertical) 
axis correspond to position (momentum) basis. Ideally, these 
Wigner functions should be nonzero only on two vertical strips 
where they take values which are ±1/8. Experimental results 
show small deviations from these values (with a maximum 
error of 15%. 



VII. CONCLUSION 

In this paper we first reviewed the formalism to rep- 
resent in phase space the states and the evolution of a 
quantum system with a finite dimensional Hilbert space. 
Then, we applied this formalism to study the phase space 
representation of quantum computers and quantum al- 
gorithms. Finally, we discussed how to perform direct 
measurements to determine the Wigner function. Our 
approach, based on the use of phase space point oper- 
ators to define the Wigner function enables us to find 
simple proofs for a variety of old and new results con- 
cerning both states and algorithms. 
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The differences between quantum and classical evolu- 
tion in phase space are clearly seen using this phase space 
representation: For the Wigner function to evolve classi- 
cally it is necessary (and sufficient) that the phase space 
point operators are transformed by the unitary evolu- 
tion operator in the same way that classical evolution 
maps phase space points. In our paper we described a 
rather interesting family of operators with this property. 
It contains all linear homcomorphisms in the torus (i.e., 
the family of all "cat" maps) which can be efficiently 
represented by simple quantum circuits (since they cor- 
respond to kicked systems where kinetic and potential 
terms act alternatively). Moreover, this family also in- 
cludes all cyclic shifts in the computational (or the mo- 
mentum) basis, reflections and the Fourier transform. At 
a deeper level all these properties are a consequence of 
the well known invariancc of the Wigner function under 
the action of the unitary representation of linear canon- 
ical transformations (that form the metaplectic group). 
Other unitary maps, that have a definite classical non- 
linear counterpart, will display for some time an almost 
classical behavior, signalled by a sharply peaked Z a p ma- 
trix. This is the semiclassical limit, that will be reached 
in the limit N — > oo, which is an important regime for in- 
teresting quantum algorithms since it corresponds to the 
limit of large number of qubits. A much wider category 
of maps can display a more drastic wave phenomenon: 
strong diffraction due to discontinuities. Consider, as the 
simplest example the phase space operation of a control- 
U gate. The controlling qubit divides phase space in two 
disjoint regions and the £/-gate acting on the remain- 
ing qubits affects these two regions differently. In this 
sense it acts as a prism or a screen in an optical system. 
Another simple example of this sort is the action of the 
operator Upj? , that docs not act on the most significant 
qbit but applies the Fourier transform to the rest (the 
diffraction effects are clearly seen in Figure 4). The im- 
plications and usefulness of these semiclassical features in 
quantum algorithms remain largely unexplored and will 
be the subject of further investigation [28]. 
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